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The energy dependence of the electron - positron hadronic cross section can be measured not only by a straight¬ 
forward energy scan, but also by means of the radiative return method. To provide extensive comparisons between 
theory and experiment a Monte Carlo event generator is an indispensable tool. We have developed such a gen¬ 
erator called PHOKHARA, which simulates e + e“ —> mesons + photon(s) processes. In this paper we present its 
latest tests and upgrades. 


1. INTRODUCTION AND PHOKHARA 
UPGRADES 

In view of the precision of the recent measure¬ 
ments of the muon anomalous magnetic moment 
= (g — 2)^/2 at BNL [1], hadronic contri¬ 
butions are crucial for the interpretation of this 
measurement, in particular for the isolation of 
the electroweak or of non-Standard Model physics 
contributions [2]. Their understanding becomes 
even more important, as it seems [3] that the 
e + e _ annihilation data are not consistent with 
the r decay data. A new measurement, which 
is under way, will challenge the theoretical pre¬ 
dictions even more. 

An important ingredient and the dominant 
source of uncertainties in the theoretical predic¬ 
tion for the muon anomalous magnetic moment 
is the hadronic vacuum polarization [4], It is 
in turn related via dispersion relations to the 
cross section for electron-positron annihilation 
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into hadrons <Jhad = u(e + e~ —> hadrons) and 
in some cases via isospin symmetry to the decay 
width r —» v T + hadrons. This quantity plays an 
important role also in the evolution of the elec¬ 
tromagnetic coupling cxqed from the Thompson 
limit to high energies [4,5]. The interpretation of 
improved measurements at high energy colliders 
such as LEP, Tevatron, the LHC or TESLA there¬ 
fore depends significantly on the precise knowl¬ 
edge of cr had . 

The feasibility of using tagged photon events at 
high luminosity electron-positron storage rings, 
such as the ^-factory, DA$NE, CLEO-C or B- 
factories, to measure <Jhad over a wide range of 
energies has been proposed and studied in detail 
in [6-8] (see also [9,10]). The machine is operating 
at a fixed energy of the e + e _ centre-of-mass sys¬ 
tem (cms) and the initial state radiation (ISR) is 
used to reduce the invariant mass of the hadronic 
system. 

The radiation of photons from the hadronic 
system (the final state radiation, FSR) should 
be considered as a background and can be sup¬ 
pressed by choosing suitable kinematical cuts, 
or controlled by the simulation, once a suitable 
model for this amplitude has been adopted. One 
finds that selecting events with the tagged pho¬ 
tons close to the beam axis and well separated 
from the hadrons indeed reduces FSR drastically. 
As demonstrated in Fig. 1, the FSR contribution 
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to the total cross section can be easily reduced to 
the 1% level. The model dependence of the FSR 
in the case of the 7 r + 7 r~ final state can be con¬ 
trolled by the same experiment through studies 
of the forward-backward asymmetry of the an¬ 
gular charged pion distribution. The asymmetry 
comes from FSR ISR interference and integrates 
to zero for ’charge blind’ configurations. As a re¬ 
sult it does not contributes to rates in Fig. 1. 
It can be used, however, to calibrate the FSR 
amplitude and more detailed tests of its model 
dependence. 
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Figure 1. The role of the cuts in the suppression 
of the FSR contributions to the cross section. Re¬ 
sults from the PHOKHARA generator. No cuts 
(upper curves) and suitable cuts applied (lower 
curves). 


When running at higher energies, the FSR is 
suppressed with respect to the ISR by the dif¬ 
ferent behaviour of various propagators and form 
factors relevant to the problem. In practice it 
means that no special angular cuts are needed to 
suppress the FSR contribution when running at 
high energies. 

The suppression of the FSR overcomes the 
problem of its model dependence, which must be 


taken into account in a completely inclusive mea¬ 
surement [ 11 ]. 

Preliminary experimental results using this 
method have been presented recently by the 
KLOE collaboration at DA$NE [12-15]. Large 
event rates were also observed by the BaBar col¬ 
laboration [16]. 

In the first version of the newly developed 
Monte Carlo program PHOKHARA [17] we have 
considered the full next-to-leading order (NLO) 
QED corrections to the ISR in the annihilation 
process e + e _ —> 7 + hadrons, for the case where 
the photon is observed under a non-vanishing an¬ 
gle relative to the beam direction. The virtual 
and soft photon corrections were presented in [18] 
and the contribution of the emission of a second 
hard photon in [17]. The final hadronic state was 
limited to the 7r + 7r _ , with the hadronic current 
modeled as in [19], and the final state emission 
was not included. The program allowed also for 
the generation of /r + /i _ 7 ( 7 ) final states, again 
limited to the emission of photon(s) from the ini¬ 
tial leptons. 

Radiative corrections proportional to (am] ) 
relevant to configurations with photons emitted 
at very low (~ m e /yfs ) angles relative to the 
beam direction were calculated in [ 20 ] and are 
included in the new version of PHOKHARA [21]. 
The leading order corrections proportional to m 2 
are typically of the order of a few per cent [ 22 ], 
while the non-leading ones are of order of 0 . 1 %, 
as seen from Fig.2. They will be important when 
the precision of the measurement will be below 
1%. Their effect depends on Q 2 and thus affects 
the Q 2 distribution from which the hadronic cross 
section is extracted. 

Another new feature of the PHOKHARA 
event generator is the inclusion of the four-pion 
hadronic final states (27 t + 27t~ and 27r°7r + 7r _ ). 
The description of the hadronic current in that 
case is based on the paper [23], with changes de¬ 
scribed in [ 8 ]. The comparison with the Monte 
Carlo, which simulates the same process at lead¬ 
ing order [ 8 ] and includes additional collinear 
radiation through structure function (SF) tech¬ 
niques, shows typical difference of order of 1 % as 
seen in Fig. 3 (a similar behaviour can be ob¬ 
served for 27r°7r + 7r _ final state). The non-leading 
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e + e —> 7r + 7r 7(7) 



Figure 2. The relative contributions of the non¬ 
leading mass corrections to the differential cross 
section at yfs = 1 GeV, 4 GeV and 10 GeV. 


contributions to the cross section of the four-pion 
final states are of the expected size and of the 
same order as for the two-pion final state [17]. 

The program now includes also the contribu¬ 
tions from the final state emitted photons to¬ 
gether with the ISR -FSR interference calculated 
at the lowest order for 7r + 7r~ and final 

states, while for the four-pion final states the FSR 
contribution is not taken into account. 

2. TESTS OF PHOKHARA 

An obvious and one of the most important 
tasks in the construction of a Monte Carlo event 
generator is to demonstrate that its technical ac¬ 
curacy is much better than the desired physi¬ 
cal accuracy. The tests that were performed for 
the previous version of PHOKHARA [17] are still 
valid, within the limitations of this version, which 
was applicable for non-vanishing photon angles. 

In the first step we demonstrate the indepen¬ 
dence of total and differential cross sections of 
the separation parameter e (called w in [17]) be¬ 
tween soft and hard photon regions. The soft pho¬ 
ton contribution is calculated analytically, while 
the additional hard photon is treated via Monte 
Carlo simulation. The parameter that specifies 
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Figure 3. The relative non-leading contributions 
to the differential cross section at yfs = 4 GeV. 
NLO - full next-to-leading result, LL - leading 
logarithmic approximation. 


the separation between the two regions of the 
phase space has to be kept small enough to val¬ 
idate the soft photon approximation and large 
enough to avoid negative weights. We performed 
the tests for a 7r + 7T _ hadronic final state in [17], 
while for one of the four-pion modes the results 
are collected in Figs. 4 and 5. From Fig. 4 
it is clear that the choice e = 10” 3 is still too 
big, whereas Fig. 5 demonstrates the stability 
of the results between e = 10“ 4 and e = 10 -5 . 
This also proves that the Monte Carlo integra¬ 
tion works well in the soft photon region. The 
Monte Carlo integration of the part of the pro¬ 
gram that simulates one hard large-angle photon 
emission was tested in [17] against a Gauss nu¬ 
merical integration. As shown in Fig. 6 a techni¬ 
cal precision of the program at the level of 10~ 4 
was demonstrated. Analytical results exist for the 
differential (in Q 2 ) cross section integrated over 
the whole angular range of the photon(s) for both 
one and two emitted photons [24]. The compar¬ 
ison with PHOKHARA can be found in Fig. 7, 
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Figure 4. The relative difference of the differen¬ 
tial cross sections for two different values of the 
separation parameter e. 


demonstrating again an excellent technical pre¬ 
cision also for the two-photon final state. The 
results presented in Fig. 7 refer to the sum of vir¬ 
tual and hard corrections to the e + e~ —> 7r + 7r - 7 
cross section, while more detailed tests can be 
found in [21]. 

3. CONCLUSIONS 

The PHOKHARA Monte Carlo event genera¬ 
tor was upgraded, allowing for simulation in the 
small photon angles region. Besides the 7r + 7r~ 
hadronic final state, its present version includes 
also 27r + 27r _ and 27r°7r + 7T _ final states. For the 
7r + 7T - and g + g~ final states, the FSR photonic 
contributions were implemented at the lowest or¬ 
der, including ISR-FSR interference. Further up¬ 
grades are in progress. 
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